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Abstract 

C*"- ' The quantum phase transition of the system of TV radiators in case of 

^3 I two-photon exchange interaction with cavity electromagnetic field is con- 

OP . sidered. It is shown that for this system the atom-field exchange integral 

increases with the increase of temperature. This effect generates pecu- 
liarities in phase transition picture. The process of nonlinear exchange 
between the quasi-spins leads to the increase of critical temperature as 
compared to traditional second order phase transition of quasi-spin sys- 
tem. 
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The cooperative emission phenomenon for dipole-forbidden transitions of in- 
verted system of radiators can be observed in the processes of two-photon spon- 
taneous emission in free space [1] or cavity [2] . It is important from the physical 
point of view to study the cooperative phenomena in a larger aspect of statistical 
physics. A great interest in the investigation of statistical and thermodynami- 
cal properties of quantum systems represents the existence of phase transition 
such as the order-disorder transition [3]. For example a radiative system like 
Dicke model [4] exhibits phase transition from the phase where the radiators 



emit non-correlated to the phase where cooperativity between the radiators is 
estabhshed in such a way that the intensity of emission becomes proportional to 
the square number of radiators, named as super-radiance. In recent years, the 
study of the phenomena which appear in physical systems when the electronic 
or atomic subsystems are in non-linear interaction with the thermostat (large 
subsystem) is of particular interest. For example, the non-linear mechanism of 
superconductivity [5, 6], bipolaron effects [7], two-photon emission [1,2,8] and 
other non-linear phenomena. Very recently, authors in [6] obtained some in- 
teresting results for the feature of phase transition for superconductivity in the 
case when a non- linear exchange interaction of the electronic subsystem with the 
phonon subsystem exists. Therefore this paper is devoted to the non-traditional 
behavior of order parameter for two-phonon superconductivity thus stating that 
the enhancing of the critical temperature becomes possible in comparison with 
usual BCS model. So we concluded that the existence of dependence of two- 
phonon exchange integral on temperature by means of the average number of 
phonons changes essentially the feature of phase transition and consequently the 
value of critical temperature with possibility to enhance this. The peculiarities 
of two-photon exchanges between two subsystems was discussed in papers [8, 
9]. In paper [9] it was shown that in special conditions the two-photon quan- 
tum oscillator can turn from incoherent to coherent stimulated emission near 
the threshold in the same fashion as the one-photon laser. Thus were obtained 
special ignition and evolution processes (similar to phase transition, see [10]) 
of two-photon laser from the incoherent to stable coherent emission in the case 
when the two-photon losses in multi-mode cavity prevail over the traditional 
one-photon losses in the resonator mirrors. 

In this paper we are interested to study the cooperative phase-transition of 
the system of A^ radiators that interact with the cavity thermostat via the two- 
quantum exchange processes. The two-photon process is possible when the one- 
photon interaction with the thermostat is forbidden. Such a possibility can be 
realized in microcavities with Rydberg atoms similarly to the experiments with 
two-photon maser effect described in [11]. For a large number of atoms in the 
cavity the phase transition through the two-photon interaction with the cavity 
modes can be possible. In this case the anomalous behavior of temperature 
dependence of the super-radiant order parameter (e.g. atomic polarization) 
is investigated. As the two-photon exchange integral between the radiators 
strongly depends on the temperature via the average number of quanta similar 
to the two-phonon superconductivity [6], so that the increasing of the order 
parameter with temperature in quasi-spin systems is observed. This effect leads 
to the increase of value of the critical temperature as compared to the traditional 
phase transition of quasi-spin system in which the exchange integral between 
the spins does not depend on the temperature. 

Let us consider an ensemble of radiators in the micro-cavity that enters into 
two-photon resonance with cavity electromagnetic field. In bad cavity limits 
when the two-photon spontaneous emission rate is less than the losses of cav- 
ity we can adiabatically eliminate the EMF operators and obtain the master- 
equation for the atomic system discussed in [12]. 
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Here the following parameters are defined: a;2iis the atom transition frequency, 

n is the average number of emitted photons in the cavity. The parameter x = 

hF' (ui '^-2uj'^)'^+iT'^ describes the cooperative atom-atom interaction through 

the cavity EMF vacuum and 7 — -jpr- 7 -2'^ X^+at'^ ^^ ^ spontaneous emission 

rate of single atom in the cavity, where A{k) is considered for the scheme of 
transitions taken into account in [12] in the form of 
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Here dijp is the moment of dipolar transition between the virtual states \a) and 
states \(3) (/3 = 1,2); g^ = yJ^jrfiUJkjV e\ is atom-cavity interaction constant, 
where e^ is polarization and V is volume of micro-cavity. In Eq.(l) are consid- 
ered the following cooperative quasi-spin operators R^ , R~ and Rz are defined 
like in [13] which are obeying the commutation relations: [i?^,i?^] ~ 2Rz and 
[Rz,R^] ^±R^. 

From the above master equation it is observed that for a larger value of the 
detuning 6 = 2a;fc — W21 comparatively to the cavity damping rate F, the value of 
two-photon decay rate, 7, is less than the absolute value of exchange integral, x- 
In this case x/l ^ 1 ^-nd we resume that the Hamiltonian part which describes 
the two-photon exchange interaction between the radiators predominates the 
losses part in the master equation (1). The effective Hamiltonian that describes 
this interaction can be expressed as 

H''" ^LjRz- XR+R^, (2) 

where the following notations are introduced : lu — h (uj2i — 2n'^x) ^n.d ^ — 
hx{i + 2n). Here A represents the energetic exchange integral between the 
radiators and it is observed that this parameter depends on the average num- 
ber of emitted photons in the cavity. We emphasize that in comparison with 
phase transition of traditional quasi-spin system like Dicke model in our case 
the exchange integral depends on temperature through the number of photons, 
n = (exp [Tiwfc/^] — 1) , — nT, where k is Boltzmann constant. Using the 
conservation of Bloch- vector at low temperature R^ — R^ — Rz + R'^R~ (where 
i? is a constant), we can express the Hamiltonian (2) through the operator Rz 

W^ = TuRz-X{R + Rz){R-Rz), (3) 

where w = uj — X = ft [uj2i ^ X ^ 2xn (1 + n)]- Thus from the condition of 
existence of the minimum value of the Hamiltonian on the operator i?^, the 
value Rz can be obtained in which the system relaxes 



and therefore the value of parameter A must be positive from the condition of 
existence of absolute minimum value of the Hamiltonian. This is possible when 
the detuning 6 takes positive values. As the ratio zu/X takes arbitrary values 
the restriction of values for parameter Rz can be obtained taking into account 
the condition of conservation of Bloch-vector, where -R^ax — j(j + !)■ From 
this expression follows that Rz must take the values between — j and j, i.e. 
—j < (Rz) < j ■ Here j = N/2 represents the maximum value of Bloch-vector, 
A'^ is the number of atoms in the cavity. Thus, for the case of a large number of 
radiators from the inequality of Rz and Exp. (4) it follows that the two-photon 
exchange integral must satisfy the inequality — 2jA < w < 2jX. The fact that 
A > is very important, because in this case the last term in the effective 
Hamiltonian (2) will imply the tendency to the order phase of the system of 
radiators. Thus, for the case when A < the phase transition is not possible. 

Let us consider that R+R^ = J2 J2k.i '^k'^T = (^z + J) +J2 J2k=ii ^t^i ■• 
where a^ are the Pauli spin operators. Here the first term describes the inter- 
action of quasi-spins with self radiation field and is similar to traditional Lamb 
shift. The last term describes the interaction of the k quasi-spin with the field 
of other I quasi-spins, while k ^ I. For large value of j this term can be ap- 
proximated in the thermodynamic limit (when lim (N/V) ~ const, while cavity 
volume is greater than the volume per atom in the cavity) in the form of 

N N N 

1=1 k=l,k=il m=l 

SO that the quadratic Hamiltonian (2) can be approximated in the manner of 
Bogolyubov method [14] by a linear model Hamiltonian. 

In the mean field approximation theory the atoms in the system can be 
regarded as identical radiators so that the spin operators ct^ do not depend on 
the atomic label. Thus the Hamiltonian (2) according to expression (5) takes 
the following form 

HA = NmJz-N\{CJ+ + C*J-)+NXCC*, (6) 

where the number-operator C = (J^) plays the role of order parameter of 
the system of radiators; here new operators J^ = i?^/-\/iV and Jz = Rz/N 
are defined which obey the commutation relations similar to the operators for 
radiators. 

Further we would like to describe the feature of phase transition of quasi-spin 
system from uncorrelated to correlated state where the cooperativity between 
the radiators is established as temperature dependence. The expression for the 
order parameter, C, can be found from the condition of absolute minimum 
of free energy since we take into account the thermodynamic analysis of the 
system. As was shown in the past in many papers related to the subject of 



phase transition for Dicke model [4] the free energy for the system described 
by Hamiltonian (2) coincides in the thermodynamic hmit with the free energy 
calculated by model Hamiltonian (6) choosing the values of parameters C and 
C* from the condition of absolute minimum of free energy. Denoting the point 
of minimum by C we have lim Fn[H'^^] = Fao[H a{C)] or in other words 



F[Ha{C)] = -(^log(Trexp[-i/^/0]) 
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Thus on condition that {dF/dC)(j^(j = 0, the value of order parameter, C, 
can be found by solving the following transcendental equation 
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The critical temperature is calculated from the consideration that for 9 



the order parameter 
solving the equation 
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and the value of this temperature is found by 
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where w and A from the above definitions arc functions of temperature and 
thus will depend on the critical temperature. Hence we see that from Eq. (9) 
it is not possible to find the explicit expression for critical temperature but the 
equation can be solved numerically. As the function hyperbolic tangent takes 
the values between —1 and 1 one can find the relations between x and W21 which 
satisfy 2fi^j. < uj2i/x < 2(1 + ncr)^ and thus the critical temperature is in strong 
connection with values of parameters x and W21. 

Exp. (8) clearly defines the important role of the two-photon exchange inte- 
gral, A, under the temperature dependence of order parameter, C. The intrinsic 
temperature dependence of the physical quantities A and w strongly influences 
the behavior of phase transition and new features are manifested in comparison 
with traditional phase transition in quasi-spin systems when exchange integral is 
constant. This peculiarity is observed from Fig.(l) in which we compare the tra- 
ditional behavior of second order phase transition with temperature dependence 
of the order parameter discussed in our model, see Exp. (8). It is an interesting 
fact that, in the case of our model the temperature evolution of the order param- 
eter in dependence with the value of the parameter x/a'21 differs substantially 
from the evolution of the order parameter calculated for traditional quasi-spin 
model which is obtained from Exp. (8) omitting the term which is proportional 
to the mean number of photons, h. Let us consider, for example, the case when 



Xl^i\ = 0.5 for which the phase transition in the usual quasi-spin model is 
absent. As follows from the numerical results of Eq.(8), the phase transition 
is possible in this case, plotted in Fig. (la). From this figure we observe that 
the order parameter appears from vacuum fluctuations of cavity field which in- 
creases with the increase of the temperature. Achieving the maximum value the 
order parameter decreases to zero value in critical point computed from Exp. (9). 
The small enhancing of value of parameter, 0.5 < xl^2i < 0.6, leads to the ap- 
pearance of phase transition in traditional quasi-spin model which exhibits a 
critical temperature to a lesser value in comparison with phase transition in 
case of proposed model (see Fig. lb, c). Further enhancing of value of x/cj2i 
will evidence inessential differences between the compared models and is of no 
interest for numerical and graphical interpretation. 

Another interesting result refers to the mean value of atomic population, 
{Rz), which can be calculated by the standard definition, (Rz) = Tr{Rzp}, 
or using the expression of free energy the following relation dF[HA]/dCj = 
— A^Tr {. . . Jz} /Tr {...} = N{Jz) = (-Rz) is observed, consequently obtaining 
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This value can be compared with the value calculated by Exp. (4) in order to 
appreciate the differences between the results obtained using the exact effective 
Hamiltonian (2) and the asymptotic Hamiltonian (6). 

From Fig. 2 wc see that values of atomic population calculated in two dif- 
ferent manners coincide in the domain of temperatures close to the critical 
temperature. This result confirms the accuracy of above calculations where the 
asymptotic Hamiltonian (6) was considered. 

In this paper we discussed the contribution of two-photon exchange mech- 
anism between the quasi-spins to the behavior of phase transition of such sys- 
tem. Therefore obtaining an intrinsic temperature dependence of two-photon 
exchange integral, the radiator system will exhibit a different fashion of phase 
transition in comparison with traditional second order phase transition of spin 
systems. The anomalous temperature dependence of order parameter in the 
case of proposed model is influenced by the appearance in the expressions of 
energy, w, and exchange integral. A, of additional terms which are proportional 
to the mean number of photons. From the numerical simulation wc obtained 
that the critical temperature can be enhanced taking into account the non- 
linear exchange process between the radiators in cavity. This increase of critical 
temperature depends on the ratio between the parameter of atom-atom inter- 
action, X, and atom transition frequency, lo2i. The calculation of average value 
of atomic population in two different ways emphasizes the agreement between 
the initial effective Hamiltonian (2) and approximation Hamiltonian (6) taken 
in thcrmodynamical limit. Hence we conclude that the model Hamiltonian (6) 
describes accurately the behavior of phase transition of quasi-spin system with 



two-photon exchange interaction where a non-traditional temperature depen- 
dence of the order parameter is manifested. 

The discussed problem may be regarded as being of academic interest due to 
the fact that such temperature dependence can be realized in more complicated 
non- linear exchanges between the spins through the thermostat. On the basis 
of two-quantum exchange mechanism between the quasi-spins, in this paper, we 
pay attention to the possibility of increasing the correlation with temperature 
in the system. 
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Figure Captions 



Fig.l The temperature dependence of order parameter for different values 
of ratio xl^2i'- a) xl^2i = 0.5; b) xl^2i — 0.51 and c) xl^2i — 0.6. Here the 
temperature is represented in relative units T/T™^^, where T™'^^ corresponds to 
the critical temperature in phase transition of proposed model. The dotted line 
corresponds to traditional phase transition of a quasi-spin system with constant 
exchange integral, A (for example Dicke model) and full line evidences the phase 
transition in the case of the proposed model. 

Fig. 2 The temperature dependencies of the average value of atomic popu- 
lation calculated using Exp. (4) - dotted hue, and by Exp. (10) - sohd hue, for 
case x/^21 ~ 0.6. 
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